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Abstract 

To study transport properties of complex networks, we analyze the equivalent conductance G 
between two arbitrarily chosen nodes of random scale-free networks with degree distribution P{k) ~ 
in which each link has the same unit resistance. We predict a broad range of values of G, with 
a power-law tail distribution <J>sf(G) ~ G~ 9G , where go = 2 A — 1, and confirm our predictions 
by simulations. The power-law tail in ^sfIC) leads to large values of G, thereby significantly 
improving the transport in scale-free networks, compared to Erdos-Renyi random graphs where the 
tail of the conductivity distribution decays exponentially. Based on a simple physical "transport 
backbone" picture we show that the conductances are well approximated by ck^ks / '(k^ + fcs) for 
any pair of nodes A and B with degrees k^ and fee. Thus, a single parameter c characterizes 
transport on scale-free networks. 
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Recent research on the topic of complex networks is leading to a better understanding 
of many real world social, technological, and natural systems ranging from the World Wide 
Web and the Internet to cellular networks and sexual-partner networks jj]. One type of 
network topology that appears in many real world systems is the scale- free network j^J, 
characterized by a scale-free degree distribution: 

P{k) ~ k~\ k min <k< A; max , (1) 

where k, the degree, is the number of links attached to a node. The distribution has two 
cutoff values for k: k min , which represents the minimum allowed value of k on the graph 



(&min = 2 here), anc 



k max = fcminiV 1 /^ which is the typical maximum degree of a network 



with N nodes The scale-free feature allows a network to have some nodes with a 

large number of links ("hubs"), unlike the case for the classic Erdos-Renyi model of random 

nn 

networks |a,lo|. 

Here we show that for scale- free networks with A > 2, transport properties characterized 
by conductance display a power-law tail distribution that is related to the degree distribution 
P{k). We find that this power-law tail represents pairs of nodes of high degree which 
have high conductance. Thus, transport in scale-free networks is better than in Erdos- 
Renyi random networks. Also, we present a simple physical picture of transport in complex 
networks and test it with our data. 

The classic random graphs of Erdos and Renyi P, |(| have a Poisson degree distribution, in 
contrast to the power-law distribution of the scale-free case. Due to the exponential decay of 
the degree distribution, the Erdos-Renyi networks lack hubs and their properties, including 
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transport, are controlled solely by the average degree 

Most of the work done so far regarding complex networks has concentrated on static 
topological properties or on models for their growth jl|, S, Isl E|. Transport features have 
not been extensively studied with the exception of random walks on complex networks 
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Il2ll. despite the fact that transport properties contain information about network 



function |l3j . Here, we study the electrical conductance G between two nodes A and B of 
Erdos-Renyi and scale-free networks when a potential difference is imposed between them. 
We assume that all the links have equal resistances of unit value jl^j ]. 

To construct an Erdos-Renyi network, we begin with a fully connected graph, and ran- 
domly remove \ — kj (N — 1) out of the N(N — l)/2 links between the N nodes. To generate 
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a scale-free network with N nodes, we use the Molloy-Reed algorithm Q|, which allows 
for the construction of random networks with arbitrary degree distribution. We generate hi 
copies of each node i, where the probability of having hi satisfies Eq. (JTJ). These copies of 
the nodes are then randomly paired in order to construct the network, making sure that two 
)reviously-linked nodes are not connected again, and also excluding links of a node to itself 

The conductance G of the network between two nodes A and B is calculated using the 
Kirchhoff method where entering and exiting potentials are fixed to = 1 and Vb = 0. 
We solve a set of linear equations to determine the potentials Vi of all nodes of the network. 
Finally, the total current I = G entering at node A and exiting at node B is computed by 
adding the outgoing currents from A to its nearest neighbors through Ylji^A — Vj), where 
j runs over the neighbors of A. 

First, we analyze the pdf that two nodes on the network have conductance 

between G and G + dG. To this end, we introduce the cumulative distribution F(G) = 
<&(G')dG', shown in Fig. QJa) for the Erdos-Renyi and scale-free (A = 2.5 and A = 3.3, 
with k min = 2) cases. We use the notation $sf(C) and Fsf(G) for scale-free, and <3>er(G0 
and .Fer(G) for Erdos-Renyi. The function Fsp(G) for both A = 2.5 and 3.3 exhibits a tail 
region well fit by the power law 

F SF (G) ~ G~( ga ~ l \ (2) 

and the exponent (go — 1) increases with A. In contrast, F ER (G) decreases exponentially 
with G. 

Increasing N does not significantly change Fsf{G) (Fig. E^b)) except for an increase in 
the upper cutoff G max , where G ma x is the typical maximum conductance, corresponding to 
the value of G at which $sf(C) crosses over from a power law to a faster decay. We observe 
no change of the exponent go with N. The increase o f GL ax with N implies that the average 
conductance G = J GQ(G)dG also increases slightly [18]. 

We next study the origin of the large values of G in scale-free networks and obtain 
an analytical relation between A and go- Larger values of G require the presence of many 
parallel paths, which we hypothesize arise from the high degree nodes. Thus, we expect that 
if either of the degrees h^ or hs of the entering and exiting nodes is small, the conductance 
G between A and B is small since there are at most k different parallel branches coming 
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out of a node with degree k. Thus, a small value of k implies a small number of possible 
parallel branches, and therefore a small value of G. To observe large G values, it is therefore 
necessary that both k A and k B be large. 

We test this hypothesis by large scale computer simulations of the conditional pdf 
$sf(G|A; j 4, k B ) for specific values of the entering and exiting node degrees kA and k B . Con- 
sider first the case k B <C and the effect of increasing k B , with kA fixed. We find that 
$SF(C|/c y i, k B ) is narrowly peaked (Fig.|2fa)) so that it is well characterized by G*, the value 
of G when $gp is a maximum. Further, for increasing values of k B , we find [Fig. Efb)] G* 
increases as G* ~ k B , with a = 0.96 ± 0.05 consistent with the possibility that as iV — > oo, 
a = 1 which we assume henceforth. 

For the case of k B ^k^, G* increases less fast than k B , as can be seen in Fig. |2fc) where 
we plot G*/k B against the scaled degree x = kA/k B . The collapse of G*/k B for different 
values of kA and k B indicates that G* scales as 



(3) 

The behavior of the scaling function f(x) can be interpreted using the following simplified 
"transport backbone" picture [Fig.|2fc) inset], for which the effective conductance G between 
nodes A and B satisfies 

- — — — 

(JA ^tb 

where 1/Gtb is the resistance of the "transport backbone" while 1/Ga (and 1/G B ) are the 
resistances of the set of bonds near node A (and node B) not belonging to the "transport 
backbone". It is plausible that Ga is linear in fc^, so we can write Ga = ckA- Since node B 
is equivalent to node A, we expect G B = ck B . Hence 

r= 1 _, ck A /k B 

l/ck A + l/ck B + 1/G tb B l + k A /k B + ck A /G tb } 1 } 

so the scaling function defined in Eq. (jSJ) is 

^ ~ l+x + ck A /G tb ^ l + x ^ 
The second equality follows if there are many parallel paths on the "transport backbone" 
so that 1/Gtb <^ ckA [19]- The prediction (jHJ) is plotted in Fig. Efc) and the agreement 
with the simulations supports the approximate validity of the transport backbone picture 
of conductance in complex networks. 
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The agreement of © with simulations has a striking implication: the conductance of 
a scale-free network depends on only one parameter c. Further, since the distribution of 
Fig. |2(a) is sharply peaked, a single measurement of G for any values of the degrees k A and 
ks of the entrance and exit nodes suffices to determine G*, which then determines c and 
hence through Eq. (0) the conductance for all values of kA and k B - 

Within this "transport backbone" picture, we can analytically calculate Fsf{G). Using 
Eq. Q, and the fact that $sf(G|A;a, k B ) is narrow, yields 

$sf(G) ~ J P{k B )dk B J P(k A )dk A 5 

where 8(x) is the Dirac delta function. Performing the integration of Eq. ((7J) using (jfjj), we 
obtain 

$sf(G) ~ CT SG [G<G max ], (8) 

where 

fi'G = 2A — 1. (9) 

Hence, for Fsf(G j ), we have F$f(G) ~ G~ 2A ~ 2 . To test this prediction, we perform simula- 
tions for scale-free networks and calculate the values of gc — 1 from the slope of a log-log 
plot of the cumulative distribution Fgp(G). From Fig. 0(b) we find that 

go - 1 = (1.97 ± 0.04)A - (2.01 ± 0.13). (10) 

Thus, the measured slopes are consistent with the theoretical value predicted by Eq. 



2l|. 



In summary, we find a power-law tail for the distribution of conductance for scale-free 
networks and relate the tail exponent gc to the exponent A of the distribution P(k). Our 
work is consistent with a simple physical picture of how transport takes place in complex 
networks. 

We thank the Office of Naval Research, the Israel Science Foundation, and the Israeli Center 
for Complexity Science for financial support, and L. Braunstein, S. Carmi, R. Cohen, E. 
Perlsman, G. Paul, S. Sreenivasan, T. Tanizawa, and Z. Wu for discussions. 
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FIG. 1: (a) Comparison for networks with N = 8000 nodes between the cumulative distribution 
functions for the Erdos-Renyi and the scale- free cases (with A = 2.5 and 3.3). Each curve represents 
the cumulative distribution F(G) vs. G. The simulations have at least 10 6 realizations, (b) Effect 
of system size on Fsf(G) vs. G for the case A = 2.5. The cutoff value of the maximum conductance 
Cmax progressively increases as N increases. 

a constant. Only for go < 2, the cutoff G max in the tail of < J>sf(C) controls G, making G an 
increasing unbounded function of N. 
[19] Flux starts at node A, being controlled by the conductance of the bonds in the vicinity of 
node A. This flux passes into the "transport backbone" , which is comprised of many parallel 
paths and hence has a high conductance. Finally, flux ends at node B, being controlled by 
the conductance of the bonds in the vicinity of node B. This is similar to the traffic around 
a major freeway. Most of the limitations to transport occur in getting to the freeway ("node 
A") and then after leaving it ("node B"), but flow occurs easily on the freeway ("transport 
backbone" ) . 
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dam, 1992). 
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FIG. 2: (a) Probability density function ^ SF (G\k A , k B ) vs. G for N = 8000, A = 2.5 and k A = 750 
(k A is close to the typical maximum degree fe ma x = 800 for N = 8000). (b) Most probable values 
of G* , estimated from the maxima of the distributions in Fig. [2(a), as a function of the degree k B . 
The data support a power law behavior G* ~ kg with exponent a = 0.96 ± 0.05. (c) Scaled most 
probable conductance G* /k B vs. scaled degree x = k A /k B for system size N = 8000 and A = 2.5, 
for several values of k A and k B : □ (k A = 8, 8 < ks < 750), (k A = 16, 16 < ks < 750), A 
(k A = 750, 4 < k B < 128), O (k B = 4, 4 < k A < 750), y ( k B = 256, 256 < k A < 750), and > 
(k B = 500, 4 < k A < 128). The solid line is the predicted function G* jk B = cx/(l + x) obtained 
from Eq. 0. This plot shows the rapid approach of the scaling function f(x) of Eq. (jHJ) from a 
linear behavior to the constant c (here c = 0.87 ± 0.02, horizontal dashed line). The inset shows 
a schematic of the "transport backbone" picture, where the circles labeled A and B denote the 
nodes A and B and their associated links which do not belong to the "transport backbone" . 
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FIG. 3: (a) Simulation results for the cumulative distribution Fsf(G) for A between 2.5 and 3.5, 
consistent with the power law Fgp ~ G~^ 9G ~^ (cf. Eq. ©), showing the progressive change of the 
slope gc — 1. (b) The exponent go — 1 from simulations (circles) with 2.5 < A < 4.5; shown also is 
a least square fit qg — 1 = (1-97 ± 0.04)A — (2.01 ± 0.13), consistent with the predicted expression 
<? G -l = 2A-2 [cf. Eq. ©]. 
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